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O ■ Abstract 

. The infrared behavior of ghion and ghost propagators offers a crucial test of con- 

finement scenarios in Yang-Mills theories. A nonperturbative study of these propa- 
^ ' gators from first principles is possible in lattice simulations, but one must consider 

\^ • significantly large lattice sizes in order to approach the infrared limit. We propose 

constraints based on general properties of the propagators to gain control over the 
extrapolation of data to the infinite- volume limit. These bounds also provide a way to 
relate the propagators to simpler, more intuitive quantities. We apply our analysis to 
the case of pure SU (2) gauge theory in Landau gauge, using the largest lattice sizes 
Qh| to date. Our results seem to contradict commonly accepted confinement scenarios. 

^ ■ We argue that it is not so. 
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■ 1 The Gribov-Zwanzieer Confinement Scenario 
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About thirty years ago, Gribov proposed an interesting confinement mechanism for color 
O I charges in Landau (and Coulomb) gauge [T]. His idea was based on the restriction of the 
^ ■ physical configuration space to the region f2 of transverse configurations, delimited by the 
so-called first Gribov horizon, where the smallest (non-trivial) eigenvalue of the Faddeev- 
Popov (FP) operator Ai = —D^d^ is zero. The limitation of the functional integration 
^ I to the (first Gribov) region fl was an attempt to fix the gauge completely, getting rid of 
I spurious gauge copies, known thereafter as Gribov copies. Since the ghost propagator G{p) 
is given by {p\M.~^\p) and Ai is semi-positive definite for gluon fields A & Q, one cannot 
have a singularity for G{p) at a finite momentum p. Using perturbation theory up to second 
order, Gribov wrote (for Landau gauge) the no-pole condition [H [2] 

Here Nc refers to the gauge group SU{Nc), an average over the Lorentz indices A has been 
considered and the quantity a{p) enters the ghost propagator as G{p) ~ [1 — o'{p)]~^. 
The above inequality tells us that, in the infrared (IR) limit, the Landau gluon propagator 
D'^'Kp) = {A'^{p)A'l{—p)) is less singular than (in the id case). Moreover, by using 
the above no-pole condition as a characterization of the first Gribov region, one can show 
that the tree-level gluon propagator becomes 

Dip) = 9'oPVip' + 7'), (2) 
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where the (Gribov) mass parameter 7 is fixed by the gap equation / d^q{2n)~*{q'^ + = 
4:/{3Ncg'^). Then, by using ([2]) and the gap equation, one can study the IR hmit of <j{p), 
obtaining an IR-enhanced ghost propagator G{p) oc l/p'^. As stressed by Gribov [1], this 
enhancement is an indication of a long-range effect in the theory that may explain color 
confinement. In Coulomb gauge, this enhancement of G{p) can be directly related to a 
color-Coulomb potential linearly rising with distance. 

It is interesting that a bound similar to the one above has also been obtained in Ref. [3] 
by considering a variational method applied to the FP operator Ai. This bound, known as 
the ellipsoidal bound, can be written (in the continuum) as 

A < c (3) 



2 



{27rY q 

where C depends on the dimensionality d of the space-time and on the gauge group. One 
should stress that the ellipsoid £, defined by the ellipsoidal bound, is a region of transverse 
configurations that includes the first Gribov region Q. A similar bound can also be obtained 
on the lattice [1]. Moreover, it is convenient to define on the lattice a region B, included in 
the ellipsoid S and including the first Gribov region Q, i.e. fl G Q G S. For all configurations 
belonging to B, one can prove [5J 

< 2tan(^), (4) 

where V is the lattice volume and A^^{0) = J2x ^^li^) is the gluon field at zero momen- 
tum. This quantity may be viewed as a magnetization M. Thus, in the infinite- volume limit 
one can show that \M\ is zero. By adding an external color "magnetic" field H coupled to A 
in the action and using the above inequality, one also obtains that the free energy per unit 
volume is null when V goes to infinity [5] . If if is spatially modulated, then the susceptibil- 
ity X at zero external field coincides with the usual gluon propagator D{p). The inequalities 
valid in the region B would then suggest that liirip^QxiH = 0,p) = limp^Q D{p) = [Ij. 

In order to restrict the functional integration to the first Gribov region Q, Zwanziger 
added to the usual Yang-Mills action a non-local term proportional to A4~^ [6]. This term 
clearly suppresses the probability of configurations near the boundary dfl of the region 
Q. After localizing the action, zeroth-order perturbation theory allows one to obtain again 
the propagator in Eq. ([2]). The purely imaginary poles p"^ = ±27^ of this propagator 
make it incompatible with a Kallen-Lehmann representation [6J. These singularities at an 
unphysical location suggest that gluons are indeed not physical excitations. One should 
also note [H [5] that -D(O) = / d'^xD{x) = means a gluon propagator in position space 
D{x) that is positive and negative in equal measure. This represents a maximal violation 
of refiection positivity. In Refs. |H [S] the violation of refiection positivity for the gluon 
propagator has been proposed as a confinement mechanism for gluons. 

The restriction of the functional integration to Q has also been discussed in Ref. [7]. 
The non-local term —7^?^ is added to the Yang-Mills action, where Ti is the so-called 
horizon function, containing an Ai^^ factor. At the same time, the Gribov mass 7 is 
fixed implicitly by the horizon condition (h) = {N^ — l)d, where h is the horizon function 
per unit volume. It is interesting that the horizon condition implies a ghost propagator 
enhanced in the IR limit i.e. limp^o[p'^G{p)]~^ = 0. Clearly, the enhancement of the 
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ghost propagator at p = should indicate that G{p) feels the singularity of M.^^ on dQ. 
Indeed, since the configuration space has very large dimensionality one expects that in the 
infinite- volume limit, due to entropy considerations, the Boltzmann weight be concentrated 
on dQ [7j. This implies that the smallest nonzero eigenvalue Xmin of A4 should go to zero 
in the infinite- volume limit. (This has been verified numerically in Landau gauge [9|.) 

The IR behavior of propagators and vertices in Landau gauge has also been studied 
in Ref. [10] by considering the sets of coupled Dyson-Schwinger equations (DSE) for the 
basic propagators and vertices of Yang-Mills theory (in the Ad case). By using a simple 
power counting for the solutions in the IR limit and constraints obtained using a skeleton 
expansion, the authors found a consistent solution characterized by an IR-enhanced ghost 
propagator G{p) ~ ^y. IR-finite gluon propagator D{p) ~ p2(2K-i) ^j^]-^ 

K e [1/2,3/4]. Note that -D(O) = for k > 1/2. On the other hand, the analysis carried 
out in [To] allows also for a solution with a tree-level-like ghost propagator at small momenta 
G{p) ~ and a finite nonzero gluon propagator D{0) > 0. These two consistent solutions 
have also been obtained by several studies of DSE using specific approximations [TT] . 

Recently it was shown [121 [13] that using the Gribov-Zwanziger approach, i.e. by re- 
stricting the functional integration to the Gribov region Q, one can also obtain in 3d and 
4:d a finite nonzero gluon propagator -D(O) and a tree-level-like ghost propagator in the IR 
limit. Here the dynamical mechanism is related to a suitable mass term that may be added 
to the action while preserving its renormalizability. As a consequence, one can show that 
the restriction to fl induces a soft breaking of the ERST symmetry [13]. It is interesting 
that the same approach cannot be extended to the 2c? case [II], because the new mass term 
produces IR singularities that make the restriction to Q impossible. 

An IR-enhanced Landau ghost propagator is also obtained as a consequence of the so- 
called confinement criterion of Kugo-Ojima [15]. At the same time, this criterion suggests 
that the perturbative massless pole in the transverse gluon propagator should disappear 
|16j . In this sense, an IR-suppressed gluon propagator (not necessarily vanishing) can 
be accommodated in this confinement scenario [17]. Finally, even though the Gribov- 
Zwanziger and the Kugo-Ojima confinement scenarios seem to predict similar IR behavior 
for the propagators, it is not clear how to relate the (Euclidean) cutoff at the Gribov horizon 
to the (Minkowskian) approach of Kugo-Ojima |18j . 



2 Bounds for the Gluon and the Ghost Propagators 

Recently we have introduced [19j rigorous upper and lower bounds for the gluon propagator 
at zero momentum -D(O) by considering the quantity 

with A'l^{0) defined in the previous section. Indeed, by straightforward calculations one 
finds that 

V{M{0)f < D{0) < Vd{N^ -1) {M{0f) . (6) 

Thus, if M(0) goes to zero as V""" we obtain that D{0) ^ 0, < D{0) < +oo, or 
-D(O) +00, respectively if a is larger than, equal to or smaller than 1/2. Recall that 
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M(0) should go to zero at least as V'^^'^ in the d-dimensional case [see Eq. At the 

same time, a necessary condition to find D{0) = is that M(0)^ goes to zero faster than 
l/V. We note that the above bounds apply to any gauge and that they can be immediately 
extended to the case D{p) with p 0. 

We investigated the bounds IQ for pure SU{2) gauge theory in Landau gauge consid- 
ering several lattice volumes in 2d, in 3d and in Ad with the largest lattice correspond- 
ing (respectively) to a'^V ^ (70 fm)^, a^V ^ (85 fm)^ and to a^V ^ (27 fm)'*. By us- 
ing the Ansatz B^/L^ for a^(M(0)^) we obtain u = 2.72(1) in the 2d case, implying 
D{0) = 0. A similar analysis in 3d and in 4d for the lower and the upper bounds gives 
0.4(1) GeV-2 < a^D{0) < 4(1) GeV-^ in 3d and 2.2(3) GeV-^ < a^D{0) < 29(5) GeV-^ in 
id. Recently, a study for the M SU{3) case [20] also finds a value for a very close to 1/2. 
Although the authors conclude that -D(O) = in the infinite-volume limit, one should ob- 
serve that in this case the lattice volumes considered are relatively small and the statistics 
is rather low. Thus, a more detailed analysis in the SU{3) case should be carried out in 
order to verify if the IR behavior of the gluon propagator agrees |21] or not with the SU{2) 
case. 

One can also obtain lower and upper bounds for the ghost propagator [22j. In Landau 
gauge, for any nonzero momentum p, one finds 

E \^mUa,p)\' < G{p) < , (7) 

where Xmin is the smallest nonzero eigenvalue of the FP operator A4 and ipminicLjP) is the 
corresponding eigenvector. If we assume Xmin ~ and G{p) ~ ^-2-2^ small p, we 
should find 2 + 2k, < 6, i.e. 5 > 2 is a necessary condition for the IR enhancement of 
G{p). Note that a similar analysis can be carried out [23] for any generic gauge condition 
J-'[A] = imposed on the lattice by minimizing a functional i?[f/], where U is the (gauge) 
link variable. Indeed, from the second variation of E[U] one can obtain the corresponding 
FP matrix A4 and the set of local minima of E[U] defines the Gribov region Q, where all 
eigenvalues of Ai are positive. In the infinite-volume limit, entropy favors configurations 
near dQ (where Xmin goes to zero). Thus, inequalities of the type ([7j) can tell us if one 
should expect an enhancement of the ghost propagator G{p) when the Boltzmann weight 
gets concentrated on dVl Q 

A study in the SU (2) Landau case [22] suggests that 5 > 2 in 2d, implying IR enhance- 
ment of G{p), while 5 < 2 in 4d. These results are confirmed if one considers the dressing 
function p'^G{p'^) for very large lattice volumes [22]. Indeed, the data in the 2d case can be 
fitted by ~ p"^", with k between 0.1 and 0.2. On the contrary, in 3d and in 4d the data 
are well described by a — 61og(l -|- cp"^), supporting k = 0. 

Let us note that our data for the gluon and ghost propagators are in good agreement 
with results obtained by other groups using very large lattice volumes [25]. Of course, one 
should also recall that the region Q is actually not free of Gribov copies [3], [5l [26] and that 
the configuration space should be identified with the so-called fundamental modular region 
(FMR) A. On the other hand, the restriction of the configuration space to the FMR should 
not make any difference on the numerical verification of the Gribov-Zwanziger scenario. 

^For example, in Ad Maximally Abelian gauge one sees that Xmin goes to zero at large volume but the 
ghost propagator stays finite at zero momentum [21] . 
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Indeed, as we have seen in the previous section, this scenario is based on the restriction of 
the configuration space to the region Q, which includes A. Actually, the bounds obtained 
for the gluon fields (see again Section [T]) apply to regions, such as 9 and £, that are even 
larger than the region Q. Finally, as explained in [27j, the restriction to the FMR can only 
make the ghost propagator less singular, as confirmed by recent lattice data [28] . 

3 Conclusions 

We have presented simple properties of gluon and ghost propagators that constrain (by 
upper and lower bounds) their IR behavior. For the gluon case we define a magnetization- 
like quantity, while for the ghost case we relate the propagator to Xmin of the FP matrix. 
We propose the study of these quantities, as a function of the lattice volume, in order to 
gain better control of the infinite-volume limit for the propagators in the IR regimeU 

Our data support a Landau-gauge gluon propagator that is IR finite in 3d and id. 
This result can be interpreted [19] as a consequence of "self-averaging" of a magnetization- 
like quantity, i.e. M(0) without the absolute value. In particular, one may think of -D(O) 
as a response function (susceptibility) of this magnetization. In this case it is natural 
to expect -D(O) > in the infinite- volume limit jj In the 2d case the magnetization is 
"over self- averaging" and the susceptibility is zero. These results are in agreement with the 
suppression of the IR components of the gluon field A due to the limitation of the functional 
space to the first Gribov region Q. At the same time the gluon propagator displays a clear 
violation of refiection positivity in Landau gauge [31], i.e. the confinement mechanism for 
gluons proposed in [H |3| is confirmed by lattice data. 

For the ghost propagator we find that in 3d and 4d the behavior at small momenta is 
essentially tree-level like, while in 2d this propagator seems to be clearly enhanced compared 
to the perturbative behavior As described in Section [H these results are not necessarily 
in contradiction with the Gribov-Zwanziger approach [T2 | [T3l [H] . 
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